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A theoretical framework is suggested for the calculation of 7*iV — »• A transition form factors 
using the light-cone sum rule approach. Leading-order sum rules are derived and compared 
with the existing experimental data. We find that the transition form factors in a several 
GeV region are dominated by the "soft" contributions that can be thought of as overlap 
integrals of the valence components of the hadron wave functions. The "minus" components 
of the quark fields contribute significantly to the result, which can be reinterpreted as large 
contributions of the quark orbital angular momentum. 
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1 Introduction 



The concept of form factors plays an extremely important role in the studies of the internal 
structure of composite particles. The non-trivial dependence of form factors on the momen- 
tum transfer (i.e., its deviation from the constant behavior) is usually a signal of the 
non-elementary nature of the investigated particle. In particular, the pioneering study of 
the nucleon form factors by Hofstadter and collaborators P demonstrated that the nucleons 
have a finite size of the order of a fermi. Later, it was observed that the behavior of the 
proton electromagnetic form factors, in a rather wide range of momentum transfers, is well 
described by the so-called dipole formula Gp{Q^)/Gp{Q) ^ GoiQ^) = 1/(1 +QV0.71GeV^)2, 
suggesting a simple Gp{Q'^) ~ 1/(Q^)^ power law for their large-Q^ asymptotic behavior. At 
the same time, strong evidence was accumulated that the pion electromagnetic form factor 
is well described by the p-pole fit F,r(Q^) ~ 1/(1 + Q'^ 1"^^ indicating that, in the pion case, 
the asymptotic behavior looks more like l/Q^. From the quark model point of view, the 
faster decrease of the proton form factor seems rather natural, since the proton contains more 
valence constituents than the pion. Furthermore, it was established that, if one can treat 
the hadrons at high momentum transfer as collinear beams of iV valence quarks located at 
small transverse separations and exchanging intermediate gluing particles with which they 
interact via a dimensionless coupling constant, then the spin-averaged form factor behaves 
asymptotically as l/((5^)^~^ [2]- This hard-exchange picture and the resulting dimensional 
power counting rules [21 E] can be formally extended onto other hard exclusive processes. 

After the advent of quantum chromodynamics, this hard-gluon-exchange picture was for- 
malized with the help of the QCD factorization approach to exclusive processes j3J E] 
that presents one of the highlights of perturbative QCD (pQCD). Within this approach, the 
hard gluon exchange contribution proves to be dominant for sufficiently large momentum 
transfers . An important ingredient of the asymptotic pQCD formalism for hard exclusive 
processes is the concept of hadron distribution amplitudes (DAs). They are fundamental 
nonperturbative functions describing the momentum distribution within rare parton config- 
urations when the hadron is represented by a fixed number of Fock constituents. It was 
shown that in the ^ oo limit, form factors can be written in a factorized form, as a 
convolution of distribution amplitudes related to hadrons in the initial and final state times 
a "short-distance" coefficient function that is calculable in QCD perturbation theory. The 
leading contribution corresponds to DAs with minimal possible number of constituents, e.g., 
3 for the proton and 2 for the pion. 

The essential requirement for the applicability of the pQCD approach is a high virtuality 
of the exchanged gluons and also of the quarks inside the short distance subprocess. Since 
the quarks carry only some fractions XjP, y^P' of the initial P and final P' momenta, 
the virtualities of the internal lines of the subprocess are generically given by XiyjQ"^, i.e., 
they may be essentially smaller than Q^, the nominal momentum transfer to the hadron. 
Assuming that (x) ~ 1/A^, one should expect the reduction factor of 0.1 for the proton and 
0.2 for the pion. In the pion case, this expectation was confirmed by an explicit calculation 
of the one-loop pQCD radiative corrections [7]. Absorbing the terms proportional to the (3- 
function coefficient (3o into the effective coupling constant ^^(/i^) of the hard gluon exchange, 
one indeed obtains {x){y)Q'^ as its argument. As a result, at accessible Q^, the bulk part 
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of the hard pQCD contribution comes from the regions where the "hard" virtuahties are 
much smaller than the typical hadronic scale of 1 GeV^ [HI IS . According to the pQCD 
factorization recipe, contributions from such regions should not be included into the hard 
term, which is strongly reduced after such contributions are subtracted. In practice, the 
subtraction is never made, and pQCD estimates are based on the original expressions (which 
implies, in particular, that the perturbative ~ behavior of propagators is trusted even 
if k"^ ^ 0). Despite this, in most cases pQCD results for hadronic form factors need special 
efforts to bring their magnitude close to experimental data. For example, assuming the 
"asymptotic" form ipn{x) = 6/7^x1X2 for the pion DA gives Q'^F^{Q'^) = Snf^ag ~ x 
0.44 GeV^ for the pion form factor ^ Ej that agrees with existing data only if one takes 
an uncomfortably large value ~ 1 for the "hard" gluon vertex. Switching to a wider 
Chernyak-Zhitnitsky (CZ) shape </?^^(x) = 30/^XiX2(xi - Xa)^ [III 1121 gives Q^F^^{Q^) = 
^vr/^as, which formally agrees with the data for ^ 0.4. However, at accessible more 
than 90% of this contribution comes from the region of gluon virtuahties below (500 MeV)^ 

i- 

In the nucleon case, the situation is even worse. For the asymptotic ~ X1X2X3 form of the 
leading twist three-quark distribution amplitude, the proton form factor turns out to be 
zero to leading order [T^ll4j . while the neutron form factor G^j is of opposite sign compared 
to the data. Assuming the equal sharing DA ~ 6{xi — l/3)5(x2 — l/3)6{x3 — 1/3) gives 
wrong signs both for the proton and neutron form factors JH]- Furthermore, if one takes 
the QCD sum rule estimate for the {qqq \ N) matrix element, the absolute magnitude of the 
form factors in the above examples is too small (by a factor of hundred) fH] compared to the 
data. Just like in the pion case, the magnitude of the formal pQCD result can be increased, 
and also the signs of the predicted proton and neutron magnetic form factors reversed to 
coincide with the experimental ones, by using CZ-type DAs [T71 UHl UHl 1201 having peaks in 
a region where the momentum fraction of one of the quarks is close to 1. Since the average 
fractions of the nucleon momentum carried by the two other quarks are small, this formal 
result is strongly dominated for accessible by regions of unacceptably small virtuahties. 
It was argued [23 that higher order Sudakov-type corrections squeeze the size of the valence 
quark configuration participating in the pQCD subprocess. Indeed, in the pion case, there are 
negative terms in the one-loop correction to the short-distance amplitude that can be written 
as Sudakov double logarithms in the impact parameter 6-space [22]- After resummation to 
all orders, they produce a factor like exp[— ln^((5^&^)/37r] suppressing the contribution of 
large transverse separations. These effects increase the region of the x, y fractions where the 
leading-order pQCD expressions are formally applicable, though they are not strong enough 
to visibly suppress nonperturbative regions for accessible momentum transfers, see e.g. [22] 
for discussion of the nucleon form factor case. 

As already emphasized, according to the standard philosophy of separating large- and 
small-virtuality contributions underlying the pQCD factorization formulas, the low-virtuality 
contributions of gluon-exchange diagrams should be treated as a part of the soft contribution. 
More precisely, in the case of the nucleon form factors, the hard pQCD contribution is only 
the third term of the factorization expansion. Schematically, one can envisage the expansion 
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Figure 1: Structure of QCD factorization for baryon form factors. 



of, say, the Dirac electromagnetic nucleon form factor Fi{Q'^) of the form (see Fig. ^ 

f.(g^) . Am + (^) ^ + (^)^ ^ + ■ ■ ■ (") 

where C is a constant determined by the nucleon DAs, while A{Q'^) and B{Q^) are form- 
factor-type functions generated by soft contributions. Because of their nonperturbative 
nature, it is impossible to tell precisely what is their large-Q^ behavior. On general grounds, 
one may expect that A^Q"^) and B{Q'^)/Q'^ correspond to higher powers of l/Q^ than the 
perturbative l/Q^ term. Perturbative estimates suggest that A{Q'^) , B {Q"^) / Q"^ < At 
very large Q^, one may also expect that they are further suppressed by Sudakov form factor. 
The most important feature of the factorization expansion is a numerical suppression of each 
hard gluon exchange by the as/vr factor, which is a standard perturbation theory penalty for 
each extra loop. If as ~ 0.3, the pQCD contribution to baryon form factors is suppressed by 
a factor of 100 compared to the purely soft term. Thus, one may expect that the onset of the 
perturbative regime is postponed to very large momentum transfers since the factorizable 
pQCD contribution 0(1/(5^) has to win over nonperturbative effects that are suppressed by 
extra powers of l/Q^, but do not involve small coefficients. In the light cone formalism, the 
functions like A(Q^) and B^Q"^) in the above expansion are determined by overlap integrals 
of the soft parts of hadronic wave functions corresponding to large transverse separations. 
There is a growing consensus that such "soft" contributions play the dominant role at present 
energies. Indeed, it is known for a long time that the use of QCD-motivated models for the 
wave functions allows one to obtain, without much effort, soft contributions comparable in 
size to experimentally observed values (see, e.g. [211 El 123)- A new trend [231211 is to use 
the concept of generalized parton distributions (GPDs, see [2H1 12^1 IHOl for recent extensive 
reviews on GPDs) to describe/parametrize soft contributions. The use of GPDs allows to 
easily describe existing data by soft contributions alone (the latest attempts can be found 
in Refs. [2111221 E2j)- A subtle point for these semi-phenomeno logical approaches is to avoid 
double counting of hard rescattering contributions "hidden" in the mo del- dependent hadron 
wave functions or GPD parametrizations. 

The dominant role of the soft contribution for the pion form factor at moderate mo- 
mentum transfers, up to ~ 2 — 3 GeV^, is supported by its calculation [211 123 within 
the QCD sum rule approach [22] applied to the vacuum average (0|T{?72(0)j(2;)?7i(|/)}|0) of 
three currents, with j representing the electromagnetic probe and the two others ?7*,?72 hav- 
ing quantum numbers of the initial and final hadrons, respectively. The application of the 
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method at higher faces the problem that the inclusion of nonperturbative effects due to 
vacuum condensates through the expansion over inverse powers of the Borel parameter 
interferes with the large-Q^ expansion of the form factors, producing a series of (Q^/M^)" 
type even for a decreasing function of Q^, like exp[— Q^/M^]. For the nucleon form factors, 
the usual QCD sum rule approach works only in the region of small momentum transfers 
Q"^ < \ GeV^ jSZlEH]- To extend the results to higher Q^, it was proposed jSHl to sum back 
the (Q^/M^)"' terms originating from the Taylor expansion of the same nonlocal condensate, 
using to this end simple models for it. Another approach [331 ED] is to use the so-called 
local quark-hadron duality approximation, in which the values of the duality intervals are 
postulated to be Q^-independent and are taken from the two-point QCD sum rules. The 
parameter-free results for the pion and nucleon form factors obtained in this way are in a 
rather good agreement with existing data. 

A less assumption-dependent approach allowing to calculate hadronic form factors for 
moderately large is based on light-cone sum rules (LCSR) [1111121 • The basic object of 
the LCSR approach is the matrix element (0|T{?72(0)j(2;)}|/ii) in which the currents r]2 and j 
are the same as in the usual QCD sum rules, while the initial hadron is explicitly represented 
by its state vector see a schematic representation in Fig. |21 When both the momentum 




Figure 2: Schematic structure of the light-cone sum rule for baryon form factors. 

transfer and the Borel parameter of the 772 channel are large, the asymptotics is 
governed by the operator product expansion T{rj2(S^)i{z)} ~ ^Cj(7/)Cj(0) on the light-cone 
= Q. The 2^-singularity of a particular perturbatively calculable short-distance factor 
Ci{z) is determined by the twist of the relevant composite operator whose matrix element 
(0|Oj(0)|/ii) accumulates nonperturbative information about the initial state parametrized 
in terms of a distribution amplitude. The lowest order O(a^) terms of the OPE correspond 
to purely soft contributions ordered by twists of the relevant operators. The magnitude 
and details of their -dependence are governed by the form of the corresponding DA. As 
shown by the LCSR calculation of the pion form factor [12], taking the CZ model for the 
pion DA, which enhances the hard contribution, one obtains a soft contribution that is too 
large. Hence one should take the DAs that are sufficiently narrower^ in order to describe the 
experimental magnitude of the pion form factor at accessible Q^. The LQSR expansion also 
contains terms generating the asymptotic pQCD contributions. They appear at proper order 

t Recent studies of the pion DAlllESlllHlilZlilllEZIEniEIlESEllEll show a converging consensus 
that the integral dx x^^ip^^ix) determining the size of the leading-order hard contribution has the value 
close to that corresponding to the asymptotic wave function, see Ref. jSS] for an updated compilation. 
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in tts, i.e., in the 0{as) term for the pion form factor, at the order for the nucleon form 
factors, etc. In the pion case, it was exphcitly demonstrated jHSl EZI that the contribution 
of hard rescattering is correctly reproduced in the LCSR approach as a part of the 0{as) 
correction. It should be noted that the diagrams of LCSR that contain the "hard" pQCD 
contributions also possess "soft" parts, i.e., one should perform separation of "hard" and 
"soft" terms inside each diagram. As a result, the distinction between "hard" and "soft" 
contributions appears to be scale- and scheme- dependent jSHl- During the last years there 
have been numerous applications of LCSRs to mesons, see |^ EH] for a review. Nucleon 
electromagnetic form factors were first considered in [001 IM] and the weak decay Af, — s> piui 
in [ng. 

In this paper, we incorporate light-cone sum rules to develop an approach to the cal- 
culation of transition form factors for electropro duct ion of the A-resonance. All three pos- 
sibilities for the virtual photon polarization are allowed in this case, hence the 7*A^ — > A 
transition is described by three independent form factors-l". Hence, the challenge is not only 
to fit the absolute magnitude of one of them, but also to explain the relations between 
different form factors. In particular, pQCD prediction IM] for the helicity amplitudes 
is Ai/2 ~ l/Q^ and A3/2 ~ 1/Q^ For multipole amplitudes Ml = -^{Ai/2 + 
and E2 = — A3/2/V3), the perturbative QCD approach predicts, hence, the same 

strength of the E2 and Ml transitions at asymptotically large Q^. Experimentally, their 
ratio is negative and extremely close to zero (within a few per cent) in the whole investigated 
region, i.e., up to ~4GeV^ jHSl EHl EZI IHHl 110] • To explain this phenomenon within the 
pQCD framework, it was suggested [TOj that the observed small value of the E2/M1 ratio 
is due to cancellation of the ^41/2 and A3/2 helicity amplitudes. It should be emphasized 
that there is no intrinsic reason inside perturbative QCD for such a cancellation to happen. 
Moreover, given the very strong sensitivity of the pQCD outcome on the form of the nucleon 
and A-isobar distribution amplitudes (e.g., by varying the CZ type DAs without changing 
their moments within the limits specified by QCD sum rule estimates, one can change the 
result for A1/2 by an order of magnitude and even reverse its sign^, see Refs. jBH EH 112]) 
such a fine-tuned cancellation looks like a miracle. Furthermore, since A1/2 and A^/2 are 
expected to behave differently with Q^, just a nearly constant value of E2/M1 ratio in a 
sufficiently wide region of (say, about 1 GeV^ wide) practically rules out the relevance of 
pQCD within such a region. 

On the other hand, the smallness of E2/M1 ratio is a famous prediction of the quark 
model. 40 years ago it was shown j72| that E2 is zero in the nonrelativistic SU(6) quark 
model, provided that quarks have zero orbital angular momentum. Small deviation of E2 
from zero was later explained either by D-wave admixtures [73J [73j or in terms of two-body 
exchange currents [70]. In the large Nc limit of QCD, it is possible to show [77| that the 
E2/M1 ratio has a smallness of order without making any assumptions about the 

angular momentum of the quarks. Small values for E2/M1 in the region < 4GeV^ were 

^The supply of possible choices and definitions of the three form factors offered in the literature seems 
inexhaustible. 

§In this connection, we would also like to remark that the model Ai/2(Q^) = ^i/2(0)/(l + Q^/Af)^ 
proposed and used in Ref. j7()| . with experimental negative value of ^1/2(0), is not able to reproduce the 
positive pQCD asymptotic results for ^i/2(Q^) quoted in Eqs. (5) and (6) of that paper. 



6 



obtained in the relativistic quark model [TH]. The 7*A^ A transition form factors were 
also calculated [7^] within the local quark-hadron duality approach motivated by QCD sum 
rules. It gives small values, within ±20%, both for the ratios E2/M1 and C2/M1 (C2 being 
the electric quadrupole or Coulombic transition form factor). Recent lattice calculations 
|8() | IRT j 1^ of the NA transition form factors up to 1.5 GeV^ gives small negative values for 
both these ratios. All these results provide strong evidence that the observed small value of 
E2/M1 has a purely nonperturbative origin. 

Another interesting feature of the 7*A^ A reaction is that the leading Ml magnetic 
transition form factor GMiQ"^) decreases with faster than the dipole fit (see, e.g., jHSj)- 
This was considered as a fact favoring pQCD since it usually gives for Gm(Q^) a result that 
is much smaller (by an order of magnitude) than that for the proton elastic form factor 
G^(Q2). In the large A^",; limit, assuming chiral and isospin symmetry, it was established 
|84[ l^ that the transition form factor Ga//(Q^) is expressed through the isovector component 
of the GPD E related to the elastic spin-flip nucleon form factor F2{Q^), which is also 
known to drop faster than the dipole fit. This observation was incorporated in Ref. jHH] to 
describe both F2{Q'^) and Gm{Q'^) using a model for the GPDs i?"''^(x, Q'^) with a Gaussian 
~ exp[— (1 — x)(5^/4xA^] plus a small power-law tail ansatz for the Q^-dependence. A Regge 
type ansatz E'^''^{x,Q'^) = e"''^(a;) -^g^g used in Ref. j^Hl- In both models, to get 

an accurate fit of the data, one needs to introduce a rescaling factor ~ 1.5 in the relation 
between elastic and transition GPDs. 

The goal of the present study is to set up the LCSR-based framework for the calculation 
of the 7*A^ — > A transition form factors. The light-cone sum rule formalism turns out to be 
considerably more cumbersome in this case because of the spin 3/2 of the A resonance. The 
local interpolating current with the A quantum numbers has also a nonzero projection on 
the spin 1/2 states with opposite parity, and the necessity to get rid of these contaminating 
contributions produces further complications. Still, as we show, it is possible to derive 
a general Lorentz decomposition and the twist expansion of sum rules for all three form 
factors in question. We explicitly calculate the leading as order sum rules and compare their 
consequences with existing experimental data. 

Apart from resolving several technical issues, our main finding in this work is that the 
soft contribution to the 7*A^ A form factors in the intermediate range is strongly 
affected by the valence quark configurations involving "minus" components of the quark 
fields that do not have the simple interpretation in terms of the leading-twist amplitude but 
rather correspond to the contributions of the orbital angular momentum |86| IH7j . The same 
conclusion was reached in (HU] for the nucleon electromagnetic form factor. In a more gen- 
eral context, large contributions of the orbital angular momentum can explain why helicity 
selection rules in perturbative QCD appear to be badly broken in hard exclusive processes 
at present energies. By construction, the LCSRs use nucleon distribution amplitudes of the 
leading and higher twist [HHl as the main input, and the results are very sensitive to their 
shape. Using the asymptotic distribution amplitudes we obtain a reasonable agreement with 
the data in the range 2 < < 6 GeV^ for the form factors. We believe that the accuracy 
can be improved significantly by the calculation of 0{as) corrections to the sum rules and 
especially if lattice data on the moments of higher-twist distribution amplitudes become 
available. A long-term goal of our study is to determine leading twist nucleon distribution 
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amplitudes from the combined fit to the experimental data on all of the existing form factors 
involving the nucleon. In this perspective, the present paper should be viewed as a step in 
this direction. 

The presentation is organized as follows. In Section 2 we introduce the necessary notation 
and set up the general framework. Section 3 contains the derivation of sum rules including 
higher twist corrections, which is our main result. The numerical analysis of the LCSRs 
is carried out in Section 4, together with a summary and discussion. The paper has three 
Appendices devoted to technical aspects of the calculation: In Appendix A we present a 
complete Lorentz-invariant decomposition of the correlation function. Appendix B contains 
the summary of asymptotic expressions for the nucleon distribution amplitudes, and in 
Appendix C the Belyaev-Ioffe sum rule for the A coupling constant is given and reanalyzed. 

2 General Framework 



2.1 Definition of the form factors 

The 7*A^ ^ A transition is described by the matrix element of the electromagnetic current 

ju = CuUluU + Cdd'y^d (2.1) 

between the nucleon state with momentum P and the A-isobar state with momentum P' = 
P — q. It can be written as 



(A(P') I MO) I iV(P)) = ApiP')Tp,^,N{P) 



(2.2) 



where N{P) and A^(P') are the nucleon spinor and the Rarita-Schwinger spinor for the 
A-isobar, respectively^. The decomposition of the vertex function 



G^iQ') [ - qp{P - q/2), + q{P - q/2)gp,] 



G3{Q^)[qpqu - q^QiBu] 



(2.3) 



defines three scalar form factors Gi{Q'^). As usual, = —q^. Following |Hn], one can also 
define the magnetic dipole Cm-, electric quadrupole Ge-, and Coulomb quadrupole Gc form 
factors instead of Gi, G21 G^: 



Gm{Q^) 



mp 



3{mp + itia] 



((3mA + mp){mA + mp) + Q 



rriA 



+ (m^ - ml)G2{Q^) - 2Q^Gi{Q^ 



Ge{Q' 



mp 



3{mp + tua) 



[m, 



m. 



G,{Q' 
rriA 



^We hope that denoting particle states by the same letters as the corresponding spinors will not create 
confusion. 
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GciQ') 



+ iml - m^)G2(g') - 2Q^G3{Q' 



3(mA + mp) [ 2 



(2.4) 



For a comparison with the hterature we also write down the form factors jHOl; Gt 
and the ratios Rem IHS] and Rsm (see e.g. jHH 1^21121) that are used in experimental papers 



2> 




Gm{Q 

W l^/AtA -t- I'l'P)' 



g^ + z/^ \ (niA + Tripy J ' ' ' 2mp 



_ i^2(g^) _ E,^ _ G^m 

Rem{Q ) - jj^^ -j^- (2-5) 




Ml(g2) Afi+ V 2mA ^^(g 



Note that there is a disagreement in the literature on the overall sign in the relation between 
Rsm and the ratio of the quadrupole and the magnetic form factors Gc{Q'^)/Gm{Q'^). We 
follow the definition from [H!^ . 

2.2 Choice of the correlation function 

To extract information about hadronic form factors within the light-cone sum rule approach 
we should analyze a matrix element in which one of the hadrons is represented by an inter- 
polating field with proper quantum numbers, and another is described by its explicit state 
vector, cf. Fig. |H Building the sum rule, we would need distribution amplitudes of the 
latter hadron. The nucleon DAs for all three-quark operators were introduced and studied 
in Ref . [SH] , while in the case of the A-isobar such an analysis is yet to be performed. Thus, 
in the present paper, we consider the correlation function given by the matrix element 

T,,{P,q) = I j d'ze''^^{Q\T{r^,{Q)j,{z)} \N{P)) (2.6) 

between the vacuum and a single-nucleon state \N{P)). The interpolating field for the 
A+-particle is taken in the form suggested in |M] 

r/^(0) = e'^'' [2(M'^(0)C7^rf^(0))M^(0) + (M"(0)C7X(0))rf^(0)] , (2.7) 
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where a, 6, c are color indices and C is the charge conjugation matrix. The contribution of 
A"*" to the correlation function in Eq. ()2.6p is given by 

Parametrizing the matrix element 

{0\v,mA{P',s)) = -^^A^;\P') (2.9) 

in terms of Aa, the coupling constant of the A^-particle to the current rj^, and using the 
standard spin summation formula for Rarita-Schwinger spinors 

E A-(POA?V') ^ -r + rn.) - - '-^ . (,,0, 

we write this contribution as 



(27r)2 ml - {P') 



1 2p;p^ p;7/3 - p;,!. 



3m/, 3mA 



r^,75iV(P). (2.11) 



This expression provides us with one of the starting points of our analysis. 

To extract the A-related information from the total correlator T^,y(P, g), we should take 
into account a subtlety that the current 77^ couples not only to isospin / = | spin-parity 
I"*" states, but also to isospin I = \ spin-parity | states. It makes sense to eliminate their 
contribution by a clever choice of the Lorentz structures. For a generic | resonance, call it 
A*, we define 

(0|r/,(0)|A*(P',.)) = ^(m.7,-4P;)A*(P',s) (2.12) 

with A* being the corresponding coupling and m* the mass. The A* spinor satisfies the usual 
Dirac equation {]/'' — m*) A*(P') = 0, and the matrix element of the electromagnetic current 
between the nucleon and a A* state takes the form 

(A*(P')|j.(0)|iV(P)) = A^(P') [(7.g^- ^g.)Pf ^*(g^) - la ,^q- ^\q')\ ^,N{P) . (2.13) 

The corresponding contribution to the correlation function in Eq. (j2.6p reads 

Tlt'\P.q) = 

= (^K^/^ - K^ S-iP'y t^^'^^'" ' ^^.aq'-F,^^] l,NiP) , (2.14) 

and the question is whether (and at which cost) the contributions in Eq. 1)2.111) and Eq. ()2.14|) 
can be separated. To achieve this, we have to understand the general decomposition of the 
correlation function in Eq. ()2.6j) into contributions of different Lorentz structures. To this 
end, it is convenient to go over to the twist decomposition in the infinite momentum frame. 
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2.3 Light-Cone Expansion 
2.3.1 Kinematics 

Having in mind tlie practical construction of light-cone sum rules that involve nucleon dis- 
tribution amplitudes, we define a light-like vector by the condition 

g.n = 0, ^2 = (2.15) 

and introduce the second light-like vector vector 

P^ = P^-\^^^^ P' = 0' (2-16) 

so that P — * p if the nucleon mass can be neglected, mp 0. The photon momentum then 
can be written as 

= + (2-17) 
F ■ n 

Assume for a moment that the nucleon moves in the positive direction, then and n~ are 
the only nonvanishing components of p and n, respectively. The infinite momentum frame 
can be visualized as the limit p"*" ~ Q — oo with fixed P ■ n = p ■ n ~ 1 where Q is the large 
scale in the process. Expanding the matrix element in powers of l/p"*" introduces the power 
counting in Q. In this language, twist counts the suppression in powers of . Similarly, the 
nucleon spinor N^{P,X) has to be decomposed in "large" and "small" components as 

N,{P, A) = W+ M N^iP, A) = iV+(P, A) + iV-(P, A) , (2.18) 

2p ■ n ' ' 

where we have introduced two projection operators 

A^ = J^ , A- = M_ (2.19) 
2p ■ n 2p ■ n 

that project onto the "plus" and "minus" components of the spinor. Note the useful relations 

^N{P) = mpN+{P) , flN{P) = ^^iV~(P) (2.20) 

Trip 

that follow readily from the Dirac equation ^N{P) = MN{P). Using the explicit expressions 
for N{P) it is easy to see that A+iV = ~ y/p+ while A'A^ = A^" ~ l/\/p+. 

The correlation function T^^, in Eq. ()2.6|) can be expanded in contributions with increasing 
twist. The leading twist contribution corresponds to the projection 

A+T++ (2.21) 

and contains the maximum power of the large momentum p"*". There exist three projections of 
the next-to-leading twist that are suppressed by one power of the large momentum compared 
to the leading twist, namely 

A_T++, A+T^+, A+T+^, (2.22) 
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and more contributions with loosing two powers of the momentum, etc. 

Each hght-cone projection in a general situation gives rise to several invariant Lorentz 
structures, in particular 



A4 




= {P 


n) 


'[A{Q\Pq)+ 4±B{Q\Pq)]^,N+{P), 


A_ 




= {P 


n] 


^[c{Q\Pq)+ i^V{Q\Pq)]^,N-{P), 


A4 




= {P 


n] 


[ ^^£{Q\Pq)+ i^T[Q\Pq)\^^N^[P), 


A4 




= {P 


n] 


{ i^Q{Q\Pq)^ ^^n{Q\Pq)]j,N+{P), 


A4 




= {P 


n] 


{J(g^Pg)+ ^^J{Q^Pq)]^,N+iP), 



(2.23) 

where e± is a two-component vector in the transverse plane that is orthogonal to q±: 

e±-q± = 0, {e±Y = l. (2.24) 

The invariant amplitudes A-J^ are not all independent. First of all, note that the n^qY^iiv- 
projection does not lead to new independent amplitudes because of the transversality condi- 
tion q^Tfj^y = 0. Second, the Rarita-Schwinger constraint = results in two relations 

g + J-mpC = 0, 

n-^I + mpV = 0. (2.25) 

Q 

Finally, more relations follow from the Lorentz symmetry. In order to find them, we write a 
Lorentz-invariant decomposition of the correlation function T^j, which involves 10 indepen- 
dent amplitudes, see Appendix A. Taking the necessary light-cone projections, we obtain 
two new relations that can be chosen as 

-mpV + mpB + 2n = 0. 
mpA-2{Pq)B-mpC-2£-2g = 0. (2.26) 

The remaining 6 independent invariant functions produce an overcomplete set of sum rules of 
the leading and next-to-leading twist for the three form factors of the 7*A^ — > A transition. 

2.3.2 The = |^ contributions 

Now, the contribution of the A+ to the invariant functions defined above can readily be 
calculated using Eq. (|2.1ip . The result has the following form: 



AA/(27r)^ 
{ml - P'2) 3ml 



A' 



2Gi + (2mA - mp)G2 + 2{mp - m/s)G: 
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J' 



AA/(27r)2 



AA/(27r)^ 

(m^ — P'^) 3m\mp 

AA/(27r)^ 1 
(m^ — ?,m\mp 

AA/(27r)2 1 



2mpGi + [mA(mA - mp) - Q'^]G2 + 2g2G3 



2(mp + mA)^! + (m^ + g^)G2 + 2[mA(mp - ttia) - g^JGa 



2(g^ - mAmp)Gi + {m^ - mp){m\ + g^)G'2 



(m^ — P'^) 6m- 



A 



2[mp{mp — m^) + Q (2mA + "^p)]G'i 



+ mA{mA - mp)^(mA + mp)G2 + 2Q^mA{mp - mA)G3 



AA/(2vr)^ 



(m^ — P'2) 6m^ 

2mAg'G3 
AA/(27r)2 g2 



2[g^ - 4mA + ("^A - 'mpy]Gi + mA(mp - m\)G2 



(m^ — P'2) 6mA 

AA/(27r)^ 1 
(m^ — P'2) 6mA 

AA/(27r)2 g2 



- 2Gi + {mp + mA)G2 + 2(mA - mp)G3 
2(3mA + mp)^! + (2mA(mA - mp) + g2)G'2 - 2g'G3 
2(mAmp - 3m^ - 2g2)]G'i 



(m^ — P'2) 6m^ 

+ (4m^(mp - mA) + g^(2mp - 3mA))G2 + 2Q'^{mA - 2mp)Gz 
AA/(27r)2 g2 



(m^ — P'2) 6m^ 



2(2mp + mA)Gi - (mA(mp + 3mA) + 2g^)G2 



+2(mA(mA - mp) + 2g^)G3 



2.3.3 The J 



p _ 1- 



contributions 



(2.27) 



For completeness, we present here the (unwanted) contributions to the same invariant func- 
tions A. . . of the negative parity spin- 1/2 A-resonances, cf. ()2.14p : 



B 



A* 



8A,/(27r)^ 
[ml - P'2) 

8A./(27r)^ 
(m2 - P'2)' 
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{ml — P''^)mp 



C 



A* 



V 



4A,/(2 



4A,/(27r)2 



2Q^Fi - m,F2 



7A* 



n 



A* 



{ml - P'2) L 
4A,/(27r)2 - 

2A,/(27r)2 
(m2 - P'2) 
2A./(27r^2 



g2(m, + mp)Fi + - ml)F2 
Q'^Fi + (mp - m,)F2 
Q'^m^Fi , 



A* 



(m2 - P'2) 

2A,/(27r)2 
{ml - P'2) 

2A,/(27r)2 
{ml - P'2) 



m*P2 



Q2 AQ^Fi-m,F2 



m,Pi + 4P2 



f2.28) 



Here Pi^2 = Ff-z^* {Q^) and m* = mA* is the mass of the A* resonance. Taking into account 
the relations in ()2.25|1 and (j2.2fjj) we find that there exist two independent combinations of 
the invariant functions that are free from contributions of = | resonances with arbitrary 
mass m*: 







I' 



. g2 ^A* _ g2^A* _ Q _ (2.29) 

On the other hand, using ()2.27|) one obtains very simple expressions for the A-isobar contri- 
butions to the same combinations: 

AA/(27r)2 



- + Q'B' 



jA + g2^A_g2^Z 



[m. 



P'2) 



AA/(27r)^ 



Q'G2{Q' 



[m. 



-g2[2Gi(Q^) + {rriA - mp)G2{Q% (2.30) 

'A ~ J 

This implies that one can construct light-cone sum rules for the 7*A^ A form factors Gi 
and G2 that are free from contamination by negative parity spin-1/2 resonances. For the 
form factor G3 this separation is not possible, unless one goes over to Lorentz structures of 
yet higher twist. Since the accuracy of light-cone sum rules is expected to decrease with 
twist, the gain of excluding spin-1/2 resonances in this case is probably not worth the effort. 
Note that the difference between the magnetic G m and electric G e form factors only involves 

2mp [(mA + mp)2 + Q2i 



Gm{Q') - Ge{Q' 



3mA {mA + rnp] 



'-G,{Q' 



(2.31) 
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3 Calculation of Correlation Functions 



On the other side, one can calculate the invariant functions A. ■ . for Euclidean virtuality 
(P')^ = (P— g)^ in terms of the nucleon (proton) DAs of increasing twist. The calculation can 
be simplified considerably by observing that only the isospin-one part of the electromagnetic 
current can initiate A^A-transitions. Hence, all correlation functions must be proportional to 
Cu — ed and a (simpler) calculation of the (i-quark contribution suffices to obtain the complete 
result. An additional simplification is due to the fact that, to the leading-order accuracy in 
QCD coupling, the two u-quarks in the A-current ()2.7|1 remain spectators and retain their 
position at the origin. It is, therefore, sufficient to retain only those terms in the general 
Lorentz decomposition of the three-quark matrix element between vacuum and the proton 
state, that are symmetric with respect to the interchange of the momentum fractions of the 
two M-quarks |88] : 

4(0|5^^VjO)«^^(0)<(^)|Ar(P)) = 

= (v, + ) (/'C)„^ (75A^), + V^mp (/'C)„^ {^^,N)^ + V^mp (7,C)„^ {Yl,N)^ 

+ V,ml (75A^), + V,ml {^,C)^^ {ta^^z^^.N)^ + V,ml (/C)„^ {^^,N)^ 

ri + (P'^cr^.C)^^ (7^5A^), + r,mp {z^ la ,.C) {^,N)^ 

+ %mp {a^,C)^p {(y^^'l^N)^ + %mp {P^a.^C)^^ {a^'^z,^,,N)^ 
+ %ml {z^'ia^.C)^^ (7^75A^)-, + {z^ P^ la ^,C) {^l^N)^ 

+T,ml {a,,C)^p {a''^ ^^,N)^ + %ml (^V^.C)^^ {^^'z,^,N)^ . (3.1) 

The expansion in (jH.lll should be viewed as an operator product expansion to the leading 
order in the strong coupling. Each of the functions Vj and 7^ is a function of the scalar prod- 
uct (P • z) and also depends on the deviation from the light-cone z"^ at most logarithmically. 
In addition, we take into account the 0{z'^) corrections to the leading- twist-3 structures, de- 
noted by and 7^*'^. The invariant functions Vi(Pz), . . . , %{Pz) can be expressed in terms 
of the nucleon distribution amplitudes Vi{xi), . . . , Ts{xi) with increasing twist, introduced in 
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twist-3 


twist-4 


twist-5 


twist-6 










Ti 


T2,n,T7 


T4, Ts, Tg 


n 



Table 1: Twist classification of the distribution amplitudes in ()3.2|) . 



Ref. 



Vi = , 
2T3 = Tr, 

2 {pz)% = -Ti + T5 + 2T8 
4(pz)T7 = T7-T8, 



2 (pz)V2 = V^-V2-Vs, 

4 (p2)V4 = -21^1 + ^3 + ^4 + 21^5 , 

4 (p2)2V6 = -V, + V2 + Vs + V4 + V5-Ve 

2r2 = Ti+T2-2T3, 

2 {pz)% = Ti - T2 - 2T7 , 

4 {pzf % = 2T2 - 2T3 - 2T4 + 2T5 + 2T7 + 2r8 

4 {pzf Ts = -Ti + T2 + T5 - Tg + 2T7 + 2T8 . 



(3.2) 



Each distribution amplitude F = Vi,Ti can be represented by a Fourier integral 



F{pz) 



-ix3{pz) 



F{xi, 1 - Xi - X3,X3) 



(3.3) 



where the functions F{xi) depend on the dimensionless variables xi,X2 = 1 — xi — X3,a:3, 
< Xj < 1 which correspond to the longitudinal momentum fractions carried by the quarks 
inside the nucleon. In difference to the "calligraphic" functions Vi{xi),Ti{xi) each of the 
distribution amplitudes Vi{xi),Ti{xi) has definite twist, see Table [TJ and corresponds to the 
matrix element of a (renormalized) three-quark operator with exactly light-like separations 
z"^ — > 0, see Table 2 and Appendix C in Ref. [SH] for the details. The higher-twist distribution 
amplitudes V2 (xj ),..., V6(xj) correspond to "wrong" components of the quark spinors and 
have different helicity structure compared to the leading twist amplitude. For baryons these 
"bad" components cannot all be traded for gluons as in the case of mesons [HHlinni- They are 
not all independent, but related to each other by the exact QCD equations of motion. As 
a result, to the leading conformal spin accuracy, the five functions V2(xi), . . . , Vei^Xi) involve 
only a single nonperturbative higher twist parameter. In the calculations presented below 
we use the conformal expansions of higher twist distribution amplitudes to the next-to- 
leading order (include "P-wave"). This accuracy is consistent with neglecting multiparton 
components with extra gluons (quark-antiquark pairs) that are of yet higher spin. Finally, 
the invariant functions Vf^ {Pz),T/^ (Pz) are twist-5 and can be calculated using equations 
of motion in terms of the nucleon distribution amplitudes jBU] . Explicit expressions for the 
distribution amplitudes are collected in Appendix A. 
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In the expressions given below we use the following shorthand notations: 

(•X3 fl—x'^ 



1 ^0 



F{x3) = I / dxiF{xi,l — xi — x'^, x' 



3J ' 



fX3 i>x'^ i>l—x'^ 

F{x3) = I dxW dx'^ dxiF(xi,l - xi -Xg',^') , (3.4) 



1 Jl JO 



where F = Vi,Ti. These functions result from partial integration in x^, which is done in order 
to eliminate the 1/pz factors that appear in the definition of nucleon distribution amplitudes 
flH.2|l . After this, the J d'^z integration becomes trivial. The surface terms sum up to zero. 
After a straightforward but tedious calculation we obtain the desired expansions: 







2 2 

x^rrip 



+ , p N4 -2^1 + V2 + V, + V, + V, + 2Ti -T2-T,- 2T-j - 2Ts) 
[xsP - qf 

2'7'"^7T7^ ^ 

+ ^^JJ'^y iYi -V2-V,-V^-V, + V,-T, + T2 + T,-T^ + 2Tr + 2Ts) 
- - ^) - ^3 - f . . . f . . f 4, (3.5) 

1 

SQCD = 4(e, - e.) y" dxs!^j^f^,{-V, + + V, + f, - f, - %) (3.6) 



+ (^1 - ^2 - 1^3 - ^4 - + n - Ti + T3 + r4 - Te + T7 + Tg- 

(X3P - g)^ 

1 1— X3 



^ :(Fi - ^2 - V^3) + , ,, (Ti - fa - T7 - + V^2 + V^s) 



(X3P - g)2 (X3P - q) 



„2^2 



Xgm 
(X3P - qY 



' ^ (^4- v^3 + ri-T2 + r3-r5-r7-2r8) 



/ 2x3ml 4x3m^g^ \ ~ ~ ~ ~ ~ ~ 
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^^^^^||^^^(f 2 - f 3 - f 4 + f 5 + f , + f s) ^ (3.7) 



1 ' 

rM{d) 





2 

A lyjj 4 ^ ^ ^ ^^X- CCS C^ C^ 

+ -V2-Vs-V^-V, + Ve-Ti + T2 + T,-T, + 2Tj + 2r8) 



(a;3P - q) 
rrip 

{X3P - qY ' {x^P-qf 



ml ^ 2ml[Q^ -4fA , _ T3 - T, + T, + T, + T3) k (3. 



1 1-2:3 





+ ^^^^,{-Ti + T3 + r4 - Te + T7 + r8)|, (3.9) 



1 l-xz 



^ ~ X ~ :x ^ 

- (^^p3^^^2 " ^3 " + T5 + T7 + T8)|, (3.10) 



1 1-X3 


~(x3P-g)4^^ +(^^p3^(-^l+^3+T7-Vi + V2 + V^3) 

2 3 

+ (^^^(^1 + T3 - fs - f, - - %^v,) 
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(X3P - qy 



(X3P - qY 
(X3P - g)6 



-(T2-T3-T4 + T5 + T7 + T8) L (3.11) 




2 

+ -T2-%- 2% - 2% + V^-V2- 2V, + ^4) 

377^2 X X X X X X 

+ / p N4 (^2 - Ts - T4 + T5 + Tr + Tg) 
(X3P - g)4 

+ (^^^^(^2 -T,-T, + T, + Tj + T,)y (3.12) 



rrtj 



JQCD = _ ejg2 y j ^^^^y, - 2Ti)(x,) 



(X3P - qy 



2 

^'""^ -(-efi + 3f2 + 3X5 + + efg + 5i/i - i/2 - 2i/3 - sv^ - w^) 



{xsP - qy 



5m 



+ / p N4 (^2 - T3 - T4 + T5 + T7 + T; 

(X3P - g)4 



2ximf 



(X3P - g)6 



(4Ti - 3T2 -T3-T4- 3T5 + ATq - 7Tj - TTg 



(X3P - g)6 
(xsP-g)^' 



T2 + r3 + r4-r5-r7-T8) ^, (3.13) 



1 1-X3 
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+ 



+ 



{xsP - qy 

(xsP - qY 

{xsP - qY 

x^nip 
{X3P - qY 



{xsP - qY 
(-f 1 + f 2 + f 5 - f 6 + 2X7 + 2f s) 

{V, -V2-Vs-V,-V, + V,) 



^2x^771^ 



+ 



{X3P - qY 

xsmpQ"^ 
{xsP - qY 

Sx^mpQ"^ 



(-8T1 + 2T2 + 6T3 + 6T4 + 2T5 - STg + IOT7 + lOTg) 



+ (xaP'-g)'^ ^^^ -V2-Vs-V,-V, + Ve)j. (3.14) 

The Borel transformation and the continuum subtraction are performed by using the 
following substitution rules: 



dx 



q{x) 



{q — xpy 



^ dx g{x) 
~(s-P'') 



^ dx 



xo 



X 



g{x) exp 



xQ^ xm^\ 



dx- 



g{x) 



^ dx q{x) 



(g - xPY Jo (s - P'^Y 



M. 



1 f'dx 



_ xm\ \ g{xo) e 



B Jxo 



X^ 



-so/Ml 



dx 



0{x) 



{q - xpy 



^ dx g{x) 
1 /"^ dx 



x^ 



+ 



Xn 



q{x) exp 



xM| M| y + xlm\ 



xQ^ xraj^^. 1 g{xo) e '^o/m^ 



xMl Ml 



2 Q2 + xlml 



q{xq) 



dxo Xo (Q^ + xlm\] 



2x0 (Q2 + x2mi)M| 



(3.15) 

where Mb is the Borel parameter, s = ^-^Q^ + (1 ~ x)m\ and xo is the solution of the 
corresponding quadratic equation for s = Sq: 



Xq 



(Q2 + _ ^2^)2 + 4^^g2 _ (g2 ^ _ ^2^) /(2mi) . 



(3.16) 
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The contributions ~ e~^o/Ml^ -g^^ ()3.15|) correspond to the "surface terms" arising from 
successive partial integrations to reduce the power in the denominators (g — xP)^^ = (s — 
P'^)^^(— x)^^ with > 1 to the usual dispersion representation with the denominator 
~ (s — P'^). Without continuum subtraction, i.e. in the limit sq — »■ cxo these terms vanish. 

In addition, in the hadronic representation for the same correlation functions ()2.27|1 one 
has to make the substitution 

\ .e-^l/M%_ (3 17) 



4 Light-cone Sum Rules 



Our general strategy is as follows. At the first step, we analyze the two light-cone sum 
rules that are obtained from the combinations of the correlation functions that are free from 
contributions of negative parity spin-1/2 resonances, cf. ()2.29p : 



^QCD _ ^QCD ^ g2^QCD 



B) 



JQCD ^ q2 ^QCD _ q2 ^QCD 



(2vr)^ 



(4.1) 



(mA-mp)G2(g2)]. 



Here B[ . . . ](M^) stands for the Borel transform with respect to P'^. From these sum rules 
we extract G\ and form factors and also Gm — Ge ~ Gi. At the second step, we determine 
the form factor G^ from the leading twist sum rule 



gQCD 



B) 



2mpGi + [m^{m^-mp)~Q^]G2 + 2Q^G^ . (4.2) 



The rationale for choosing this structure is that in the sum rule for A the contribution of 
interest is multiplied by a small factor ttia — mp, cf. ()2.27j) . Finally, we rewrite our results 
in terms of G^, Ge and Gc- 



4.1 Asymptotic expansion 

As an example, consider the contribution of the leading-twist nucleon distribution amplitudes 
to the LCSR for G2{Q'^), the first equation in ()4.1|) . Putting everything together, one obtains 
to this accuracy 

1 



+ (^|^(^i-^i)}+^(twist-4). (4.3) 

The main effect of the continuum subtraction and/or Borel transformation is that the inte- 
gration over the momentum fraction X3 gets restricted to a narrow interval of 1 > X3 > Xq 
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with 1 — xq = sq/Q^ <^ 1. As the result, power counting of 1 — X3 factors in the integrand 
translates into the power counting in The behavior of the nucleon distribution am- 

plitudes close to the end point X3 — > 1 is governed by conformal symmetry [HZj. To the 
next-to-leading order accuracy in the expansion over the conformal spin one obtains fZ| |HH] 

Vi{xi, X2, X3) = 120fNXiX2X3 [1 + ^(1 - 3Vi'^)(l - 3x3) + . . . ] , 

ri(xi,X2,X3) = I2O/JVX1X2X3 [1 + ^ (3^5' + 31^1'^ -1) (1-8x3) + ...]. (4.4) 

Here A" and are scale-dependent parameters that characterize the deviation of the dis- 
tribution amplitude from its asymptotic form at ^ 00. For the asymptotic distribution 
amplitude Vf = 1/3 and A" = 0; the Chernyak-Zhitnitsky (CZ) model ^7] corresponds to 
Vi = 0.23 and A" = 0.38 at a low scale of a few hundred MeV. A simple calculation shows 

that in the end-point region Vi ~ Ti ~ (1 — ^3)^ and V^i ~ Ti ~ (1 — ^3)^. 

After some algebra we obtain the contribution of leading-twist nucleon distribution am- 
plitudes to the light-cone sum rule ()4.H) for G2{Q'^) in the — > 00 limit: 

^ G'r'm = jdse-/^'^'^ J dx,V^{x„s/Q^-x„l-sm 



SOfNmp 





so 



Q 



10 



[1 + 7(31^'^ - 1)] J ds e-'l^'B + o . 



(4.5) 



The term in ~ (3Vi^ — 1) corresponds to the contribution of subleading conformal spin, 
cf. fl4.4|) . For the CZ-model this correction is very large, [1 + li^V^ — l)]c2 — —1-17. It 
turns out, however, that the contribution of the leading-twist distribution amplitudes is only 
subleading for large Q^, and the true large-Q^ asymptotics of the sum rules is determined 
by the contribution of higher-twist amplitudes, in particular those corresponding to different 
helicity structures (compared to leading twist). Using full expressions, i.e. taking into 
account the terms ~ V2 . . . Ve and ~ T2 . . . Te, we obtain 



so 



''^(y° G.(Q') = ^15/. + /d..= e-« + O (^) . (4^6) 



[2G,{Q') + {m^-mp)G2{Q^)] = ^^^0/^ / dss^e-^'^^B 



89 



In + 7Ai 



dss^e-'/''"A + 



Q 



10 



(4.7) 
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(27r)2 



so 



[75fN + 9Ai] / ds s 



s/M 



+ mp 
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34/;v + -Ai 
5 



O 



Q 



12 



(4. 



from the three sum rules in Eqs. ()4.H1 . (j4.2j) . respectively. For simplicity, only the contribu- 
tions of asymptotic distribution amplitudes are shown, see Appendix B. 

To avoid misunderstanding, note that these expressions correspond to the "soft" contri- 
bution to the form factors that does not involve any hard gluon exchanges and correspond, 
loosely speaking, to the first term in the expansion in Eq. p.lj) . In order to observe the 
true asymptotic behavior one has to calculate radiative corrections to the sum rules, cf. jSE]- 
In this way, the hard gluon exchange contribution considered in [H3] appears to be a part 
of the two-loop 0{a1) correction. The corresponding calculation goes beyond the scope of 
the present work. In spite of being enhanced asymptotically by powers of the momentum 
transfer, the radiative corrections are accompanied by increasing powers of the (small) QCD 
coupling at a scale at least that of the Borel parameter, and we expect that for moderate 
momentum transfers in the range 1 < < 10 GeV^ their contribution is numerically less 
important. 




Figure 3: The average value {x) of the momentum fraction X3 in the sum rules (|4.H) . (|4.2)) as a 
function of the momentum transfer Q^. The lower and the upper soUd curves correspond to the 
first and the second sum rule in (|4.1|) . respectively. The dashed curve corresponds to (|4.2j) . The 
value of the threshold xq is plotted by dots for comparison. 

We conclude that the leading-order LCSR calculation predicts a universal falloff of 
soft contributions to all the three form factors Gi, G21 G^, and to this accuracy G2 = 2G3. 
This translates to the asymptotic behavior of the soft terms Gm ~ l/Q^? G^ ~ 
and Gq ~ ^/Q^, for the magnetic, electric and quadrupole form factors defined in Eq. (12. 4p . 
respectively. In agreement with the common wisdom, soft contributions in the light-cone sum 
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rules arise from the integration regions where the quark interacting with the virtual photon 
carries almost all hadron momentum, alias the corresponding momentum fraction x 1. 
To illustrate this feature, we have plotted in Fig. IHlthe average value (x) of the momentum 
fraction X3 in the integrals in the sum rules ()4.H1 . ()4.2|1 . The main effect here is that the 
integration region in the momentum fraction gets restricted to a narrow interval xq < a; < 1 
where Xq is given by Eq. ()3.16|) . For asymptotically large one finds Xq ^ I — Sq/Q'^ and the 
integration region shrinks to the end-point. For realistic values of in the range 1 — 10 GeV^ 
the average (x) grows slowly and one finds very similar values for all the three sum rules in 
question. [The irregular behavior that is seen for the sum rule ()4.2|1 around ~ 2 GeV^ is 
due to accidental vanishing of the denominator involved in taking the average.] 

4.2 Numerical analysis 

The numerical results presented below are obtained using two models of the nucleon distri- 
bution amplitudes [HE]- The first model corresponds to asymptotic distribution amplitudes 
of all twists. The corresponding expressions are collected in Appendix B. The second model 
corresponds to taking into account the corrections to the asymptotic distribution amplitudes 
that are due to the next-to-leading conformal spin ("P-waves") with parameters estimated 
using QCD sum rules. For the leading twist distribution amplitude this choice corresponds 
to a simplified version of the Chernyak-Zhitnitsky (CZ) model: we have truncated the orig- 
inal CZ expressions ^Tj leaving out contributions of the next-to-next-to-leading conformal 
spin operators ( "D-waves" ) in order to simplify the calculation of nucleon mass corrections, 
see jSHl for details. The explicit expressions of the distribution amplitudes including P-wave 
corrections are rather cumbersome, they can be found in |88t l98] . 

In addition, for the numerical evaluation of the sum rules we have to specify the values 
of the A-coupling Aa, the continuum threshold sq and the Borel parameter M^. The usual 
strategy is to determine all of them from the simplest QCD sum rule involving the A reso- 
nance ^ni; see Appendix C. We obtain from this sum rule Sq — 3 GeV^ and Aa — 2.0 GeV^. 
Both numbers are somewhat lower that the ones obtained in because we use the experi- 
mental input value of the mass of the A- resonance ttt-a = 1.23 GeV instead of trying to find 
it from the sum rule. 

A suitable range of Borel parameters for the LCSR can be obtained as follows. On the one 
hand, has to be small enough in order to guarantee sufficient suppression of higher mass 
resonances and the continuum in the hadronic representation for the correlation function. 
This is the same criterium that is applied to the Belyaev-Ioffe sum rule ()C.1|) for the coupling 
Aa, hence we would want to take as close as possible to 1.5 GeV^ which is the minimum 
value at which the stability in ()C.1|) sets in. On the other hand, the Borel parameter in 
the LCSRs has to be large enough to guarantee convergence of the twist expansion in the 
QCD calculation. Note that for a fixed value of the momentum fraction x the light-cone 
expansion of the relevant correlation functions goes generically in powers of l/(g — xP)"^ 
which translates to the expansion in powers of l/(a;Af^) after the Borel transformation. 
The true expansion parameter is therefore of order ~ l/((x)M^) where (x) is the average 
value of the momentum fraction in the corresponding integrals, rather than 1/M^ itself. 
In the range 1 < < 10 GeV^ we find 0.4 < (x) < 0.8 (see Fig. Ej) so that preferred 
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Figure 4: The Borel parameter dependence of the ratio Ga/ /SGdipoic where Gdipoic = 1/(1 + 
(5^/0.71)^. The magnetic form factor Gm is defined in ()2.4() . The solid and the dashed cm'ves 
correspond to the calculation using M| = 2.5 GeV^ and M| = 1.5 GeV^, respectively. In both 
cases the asymptotic nucleon distribution amplitudes of the leading and higher twists are used. 

values of in LCSRs appear from this side to be a factor 1.2-2.5 larger compared to the 
two-point sum rule Borel parameter M^. Note that for fixed there is effectively a bias 
towards using larger values of with decreasing Q^. All in all, it appears that the interval 
1.5 < < 2.5 GeV^ presents a reasonable choice. The stability of the LCSRs in this range 
turns out to be better than 10-15%, see Fig. 0] for an example. In what follows we use the 
fixed value of the Borel parameter = 2 GeV^ in the center of this fiducial interval. 

The results for the magnetic transition form factor are shown in Fig. Fig. Inland Fig. [7| 
using three different ways to present the data (in terms of G^f^, Gm and Gt) accepted in 
the experimental papers. The difference between G^/^ and Gm is in a kinematical factor 
only, while Gt includes in addition the contribution of the electric form factor Ge which is 
negligible. We observe that the calculation with asymptotic distribution amplitudes (solid 
curves) is much closer to the data so that large deviations from the asymptotic shape of 
distribution amplitudes suggested by QCD sum rule calculations are disfavored. A similar 
conclusion was reached in Refs. [^01 1^ from the LCSR analysis of the electromagnetic 
nucleon form factors. 

It turns out that the sum rules for Gm are dominated by contributions of subleading 
twist-4. To illustrate this issue, we have plotted in Fig. |S1 separate contributions to the 
sum rule for the ratio G m / {"^G dipoie) that come from the nonperturbative matrix element 
of twist three ~ /at (dashed curves) and the two existing matrix elements of twist four: 
terms in Ai (dotted curves) and A2 (dash-dotted curves), cf. Appendix B. The full result is 
given by the sum of all three contributions and is shown by the solid curves. We see that 
contributions ~ A2 are numerically very small, the contributions ~ Ai are the dominant ones 
and remain roughly constant in the considered range, whereas the contributions ~ 
have a stronger dependence and enter with an opposite sign. QCD sum rule motivated 
corrections to the asymptotic distribution amplitudes generally tend to increase both leading- 
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Figure 5: The ratio Gff /SGdipoic where Gdipoie = 1/(1 + QVO.TI)^. The form factor Gjf is 
defined in (|2.5j) . The solid and the dashed curves correspond to the calculation using the asymptotic 
and the QCD sum rules motivated nucleon distribution amplitudes of the leading and higher twists, 
respectively. The data points are from |SS] (blue squares) jlUUj (red triangles) and |1U11 11U21 11U31 
11041 inn] (red squares). 

and higher-twist contributions, so that the cancellation becomes more pronounced. Note 
that the uncertainties in nonperturbative parameters /at, Ai, A2 presently are at the level of 
20-30%. Larger leading twist contributions ~ /at would yield a steeper falloff of the form 
factor which is favored by the data. The same observation was made in [001 ^^^^ 
of the electromagnetic form factors of the nucleon. Note that the dominant contributions 
~ Ai correspond to the operators that include a "minus" component of one of the valence 
quark fields. They do not have a simple partonic interpretation in terms of quark parton 
amplitudes at small transverse separation but rather involve the orbital angular momentum, 
see p¥)| W>7\ for the relevant formalism and discussion. 

The results for the electric form factor Ge and quadrupole form factor Gc are shown in 
Fig. Eland Fig.^1 respectively. In both cases we plot the experimentally measured quantities 
Rem and Rsm f|2.5|) that are related to the form factor ratios, normalized to the magnetic 
form factor Gm- Here, again, the asymptotic distribution amplitudes tend to give a better 
description. 

In the future, one should try to constrain the parameters in the distribution amplitudes 
by making a combined fit to the LCSR for the electromagnetic and weak form factors of 
the nucleon, including A-resonance production etc. In order to make this program fully 
quantitative one first has to calculate radiative corrections to the LCSRs, similar as this has 
been done for the pion form factor [SEl and B- meson decays [001 • The corresponding task 
goes beyond the scope of this paper. 
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Figure 6: The ratio Cm /SGdipoic where Gdipoie = 1/(1 + Q^/0.71)^. The magnetic form factor Gm 
is defined in (|2.4|) . The identification of the curves and data points is the same as in Fig. [3 



5 Summary and conclusions 

In this paper, we incorporated light-cone QCD sum rules approach to calculate the purely 
nonperturbative soft contribution to the 7*p A"*" transition. The soft contribution cor- 
responds to the so-called Feynman mechanism of the large momentum transfer, which is 
suppressed asymptotically by power (s) of compared to the perturbative contribution 
of hard rescattering. We argued from the very beginning that, due to two hard gluon ex- 
changes required in the pQCD contribution, the latter is numerically suppressed by a factor 
of 1/100 compared to the soft term. Indeed, our results for the dominant magnetic form 
factor Gm{Q'^) are rather close to the experimental data in the region above ~ 2 GeV^. 
This confirms the general wisdom that the soft (end-point) contribution is dominant at the 
experimentally accessible momentum transfers. Moreover, the inspection shows that the soft 
contribution is dominated by quark configurations involving a "minus" light-cone projections 
of one of the quark field operators, which can be reinterpreted as importance of the orbital 
angular momentum (cf. |HUl I87j ). 

In the region of low < 2 GeV^ our results for the magnetic form factor appear to 
be factor two below the data. One reason for that are most likely the so-called "bilocal 
power corrections" that correspond to long-distance propagation in the channel. Such 
terms were studied in case of the pion |106j and nucleon |SZ] form factors, and they provided 
a sizable enhancement of form factors at low momentum transfers as compared to the ex- 
trapolation of the large-Q^ results. Another possibility is that the interpolating current for 
the A-particle is not good enough and couples strongly to the excited states. This can be 
checked by applying the light-cone QCD sum rules to the correlation function (A| J?77v|0) of 
the electromagnetic current and a local current ?7jv with nucleon quantum numbers, while the 
A is left explicitly in the final state (A|. To pursue such a program, one needs a systematic 
analysis of the distribution amplitudes of the A as the first step, the task which is interesting 
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Figure 7: The ratio Gr/SGdipoie where Gdipoie = 1/(1 + Q^/0.71)^. The form factor Gt is defined 
in (|2.5)) . The identification of the curves and data points is the same as in Fig. El 

in its own. 

Our calculations are in agreement with the experimental observation that the ratios 
E2/M1 and C2/M1 are small. The charge form factor Ge appears to be very sensitive to 
the shape of the nucleon distribution amplitudes so that the experimental data can easily be 
fitted by assuming moderate corrections to the asymptotic form, of the same sign but much 
smaller as compared to the CZ model. However, both the GEiQ"^) and Gc(Q^) form factors 
in our calculations appear as small differences between expressions dominated by the much 
larger form factor GmIQ"^), a situation similar to that encountered in the approach [TU] based 
on the analysis of the 3-point correlator. For this reason, at this stage we restrict ourselves to 
a conservative statement that GeIQ"^) and Gc{Q'^) are small compared to Gm{Q'^) without 
insisting on a specific curve (or even sign) for them. 

The inclusion of QCD radiative corrections to the sum rules presents another impor- 
tant issue, and is needed to make the theoretical studies of the 7*p — > A"^ transition fully 
quantitative. Such corrections contain both purely soft contributions and also hard gluon 
exchanges, with nontrivial interplay and connections between these two types of dynamics. 
For example, for the pion form factor it was found jSH] that there exists a partial cancel- 
lation between soft contributions and hard contributions of higher twist. A big advantage 
of the light-cone sum rule technique is that it is free from double counting: soft and hard 
contributions can be separated rigorously. In the present case, such a separation becomes 
necessary starting at two-loop 0{a1) corrections to the sum rules. 

To summarize, we believe that the light-cone sum rule approach currently offers the 
best compromise between theoretical rigor and the applicability to present and planned 
experiments involving elastic and transition form factors for baryons. This approach is 
rigorous as far as the separation of hard and soft dynamics is concerned, and provides one 
with a useful tool for the study of the transition region between hard perturbative and 
soft nonperturbative QCD dynamics. One goal of such studies is to determine nucleon 
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Figure 8: Separate contributions to the light-cone sum rule results (solid curves) for the ratio 
GM/3Gdipoio where Gdipoic = 1/(1 + (5^/0.71)^. The terms in /at, Ai and A2 are shown by the 
dashed, dotted and dash-dotted curves, respectively. The results on the left panel correspond to 
the asymptotic distribution amplitudes and the ones on the right panel are obtained including a 
CZ-like model, see text. 



distribution amplitudes from the data on form factors, similar as parton distributions are 
extracted from the measured deep inelastic structure functions. Our work presents a step in 
this direction. 
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A Lorentz-invariant decomposition 

The correlation function Tn^{P,q) has to satisfy two conditions 

7%,(P, g) = , g%.(P, g) = , (A.l) 

that follow from the property of the current r/^ and electromagnetic current conservation, 
respectively. In addition, T^^^P^q) has to be proportional to the nucleon spinor. This 
suggests to define 

TV 

T^P^q) = Y,UP'^q')L^;:XP,qHN{P) (A.2) 

k=l 

where Tk{P''^, q^) are scalar functions in front of the Lorentz structures L^I^}{P' , q) that satisfy 
the conditions in Eq. themselves, i.e. YL^^Xp' ^(l)N{P) = and q''L^^J{P',q)N{P) = 

0. The complete basis contains ten different Lorentz structures which can be chosen as 

^2 = 2(g^7^^ - gf,^q^) + 7^(7^?^ - qu^) , 
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Figure 9: The ratio -Rem(<3^) = —Ge/Gm of the electric and the magnetic form factors, cf. H2.5|) . 
The solid and the dashed cm'ves correspond to the calculation using the asymptotic and the QCD 
sum rules motivated nucleon distribution amplitudes of the leading and higher twists, respectively. 
The data points are from jlU5j (red squares) and (HHI (blue squares). 

= A{q,Pl - g,,P'q) + ^,{^,P'q - P'J) , 
Lfl = P^lJ - Qy) - \mp{ljlu - quip,) + ^li^iluQ^ - Qui) + ]pf,{-iuP'q - Pli) , 

= P'qP'M - + - P'.Ki) + P"iQ,Pl - 9,.P'q) , 

^IS = - 9^^uq^ + ^7^(7^?^ - Qu^) (A.3) 

and 

= 2(g^7^ - gf,J) - {-fj-f,y - g^7;,) , 
L(Jj = 4(g^P;^ - g,,P'qi) + 2P'gg,7^ - P'q-fJ-f, - Pl^.q" , 

L^'l = P^M - Qui) + P',mp{^J - q,) + ^mp{^^^,P'q - ^^P'J) - \p'\l4lu - q.l,) 

+^{'^P'<iif^(iu - ipP'qiiu - Plq^if,) , 

= Qt^qJ - Ot^uq^i - ^q^iljlu - It^qu) , (A.4) 

= mpP;(7.^ - q.) + P;(7.g' " qj) + 2(P'gP;7. " P'.Pl^i) - \p"Mlu - q.l,) , 

where the two groups roughly correspond to contributions with even and odd number of 
gamma-matrices, respectively. 



0.7^ 
0.5 
.25 



-0 .25 
-0 . ^ 
-0.73 



30 




QiGeV 



Figure 10: The ratio Rsm ~ —Gc/Gm of the quadrupole and the magnetic form factors, cf. H2.5() . 
The solid and the dashed curves correspond to the calculation using the asymptotic and the QCD 
sum rules motivated nucleon distribution amplitudes of the leading and higher twists, respectively. 
The data points are from jlU5j (red squares) and (HHI (blue squares). 



Taking the necessary light-cone projections we obtain for the amplitudes in ()2.23|) : 

A = Q'T4, 

B = -mp T4 - 2 Tio , 
mpC = -Q^[mpTi + 2 T7 - Tg + 2 T^] , 
mpV = -2T2 + T3 + {2P'q - ml) T4 + mp Tg , 

£ = mp{P'q) T4 - Q2 Tg + {2P'q - Q^) T^, , 
= T3 + {P'q) T4 + mp Tg + mp Tio , 

H = -T2 - ^Ts - ^mp Tg , 

J = -Q2[3T2 - - {2P'q - m^) T, - ^mpTg] , 

J = + 32^^ _ + 2Tio] . (A.5) 

Note that four invariant amplitudes Ti, T5, Tg, Tg do not contribute to the next-to- leading 
twist accuracy. This implies that the invariant functions in ()2.23|) obey four relations which 
can be chosen as in Eqs. ()2.25|) and ()2.26p in the text. 
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B Nucleon Distribution Amplitudes 

To make the paper self-contained, we collect here the necessary information on the nucleon 
distribution amplitudes that enter the sum rules. Our presentation follows Ref. [HE]- 

A standard tool to study hadron distribution amplitudes is to use constraints on oper- 
ator mixing and equation of motion relations that originate from the conformal symmetry 
of the QCD Lagrangian ^7\. This approach suggests an expansion in contributions of in- 
creasing conformal spin. The leading-order contributions with the lowest conformal spin are 
usually referred to as asymptotic distribution amplitudes. To this accuracy there are only 
three nonperturbative parameters involved. One obtains [SHI asymptotic twist-3 distribution 
amplitudes: 

Vi{Xi) = 120XiX2XsfN , Ti(Xi) = 120XiX2X3/Ar , (B.l) 

asymptotic twist-4 distribution amplitudes 

V2{xi) = 2xia:2 [5(1 + X3)/Ar + 6Ai] , T2{xi) = 4x1X2X2 , 
Vsixi) = xs [5 (1 + 2x1X2 - 4x3 + 3(a;2 + xl + xj)) /jv - 6(1 - a;3)Ai] , 
Tsixi) = X3 [5 (1 + 2xiX2 + 2x3-3 [xl + xl + Xg)) /a? + (1 - a;3)A2] , 
T7(x,) = X3 [5 (1 + 2x1X2 + 2x3-3 {xl + xl + xl)) /tv - (1 - a;3)A2] , (B.2) 
asymptotic twist-5 distribution amplitudes: 

V^ixi) = [28 - 65{xl + xl) - 30x1X2 - 13x3 - 15x^] /^v - [l + (xi - X2)^ - x^] Ai , 

T4(xi) = --[28 - 65{xl + xl) + 30x1X2 - 43x3 + 15x^] /jv 
3 

+ Tl + x^ - 6x1X2 + xl + 2x3 - 3x3] A2 , 



Ts{xi) = -i [28 - 65{xl + xl) + 30xiX2 - 43x3 + 15x^] 
- ^ Tl + Xl - 6x1X2 + xl + 2x3 - 3x3] A2 , 



y5{Xi) = -^X3(28 - 65X3)/jV + 2X3(1 + X3)Ai , niXi) = ^X3(l + X3)A2 , (B.3) 

asymptotic twist-6 distribution amplitudes: 

V&ix^) = 1(5 + 3xs)fN - 1(1 - 3x3)Ai , T6(x,) = J(8 - 6x3)/jv , (B.4) 
3 5 3 

and, finally, the x^-corrections to the same accuracy: 

Vf^''\x3) = ll [(1 - ((215 - 529x3 + 427x^ - 109x^) + 41n[x3]) /n + 16(1 - X3)'Ai] , 

<^'^(^3) = f (1 - x-s)' [(76 - 39x3)/^ - 6A1] , (B.5) 
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where 



Jo 

T^^{Pz) ■= C dx3 e-*"^(^") Tf^''\x3) , (B.6) 
Jo 

cf. dUD, (unD. 

The nonperturbative parameters /at, Ai and A2 correspond to the vacuum-to-nucleon 
matrix elements of the three-quark local operators [SH] • The existing estimates come from 
QCD sum rules. At the scale = 1 GeV one gets 

fN = (5.0 ± 0.3) ■ lO-^GeV^ dH 1201 . 
Ai = -(2.7 ± 0.5) ■ lO-^GeV^ jSHl EH] , 

A2 = (5.4± 1.0) ■ lO-^GeV^ IHHl EH- (B.7) 

Calculations in this work are done also using a more sophisticated model for the distribution 
amplitudes, taking into account corrections that correspond to the contributions of operators 
with the next-to-leading conformal spin j^T] . The corresponding expressions are given in |HH1 
(see also [SHI). They include five more parameters (the first two leading twist, and the rest 
twist-4) which we choose as 

= 0.38 ±0.15 im, 
= 0.23 ± 0.03 [m , 
= 0.40 ± 0.05 [HHl EHl , 

= 0.22 ± 0.05 [Ha IHHl , 

/i" = 0.07 ± 0.05 jHl Ei • (B.8) 

For a comparison, the asymptotic distribution amplitudes correspond in this more general 
parametrization to = 0, Vf = 1/3, ff = 3/10, /| = 4/15 and = 1/10, see jHHlEH] for 
details. 

Note that there is a mismatch between the twist classification of the distribution ampli- 
tudes that implies counting powers of the large momentum p+, and the twist classification 
of the local operator matrix elements: E.g. the parameters of the twist-4 distribution am- 
plitudes depend both on the the leading twist-3 matrix element /at and the twist-4 matrix 
elements Ai, A2. Such contributions of lower- twist matrix elements to higher-twist distribu- 
tion amplitudes are well known and usually referred to as Wandzura-Wilczek contributions. 
In particular the distribution amplitudes of twist-5 and twist-6 are entirely of Wandzura- 
Wilczek type since there exist no geniune twist-5 and twist-6 operators to this order in the 
conformal expansion. 
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C QCD sum rule for A-resonance 



The coupling Aa (j2.9|l can be found from the Belyaev-Ioffe sum rule jTHj for the correlation 
function of the two 77^ currents. The sum rule reads 



2 -ml/AP 
A 



5 



+ 



4a2 



1 - e-^°/^^' 



1 

2M4 



— Ar 
72 



9m| 



where is the Borel parameter, sq the continuum threshold and 



a 
b 



rrin 



{2nY {qq) - 
2 / ^(^2 

TT 
2 



(27r) 

0.8 GeV 



0.55 GeV^ , 
= 0.47 GeV^ 



(C.l) 



(C.2) 



Using an experimental value mA = 1232 MeV we found that the sum rule (jdljl is stable in a 
broad interval of Borel parameters > 1.5 GeV^ for the choice of the continuum threshold 
So = 3.0±0.2 GeV^. The corresponding value of the coupling is equal to Aa = 2.0±0.2 GeV^. 
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